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We study the effects of fermiopliobic scalar /pseudo-scalar oblique corrections on bound state en- 
ergy levels in muonic atoms. To make the treatment sufficiently general, while including ordinary 
scalar and axion-like pseudo-scalar fields as special cases, we consider Unparticle scalar /pseudo- 
scalar operators with couplings predominantly to photons. We derive the relevant vacuum polar- 
ization functions and comment on the functional forms of the Unparticle Uehling potentials for 
various scaling dimensions in the point nucleus and finite nucleus approximations. It is estimated 
that for an infra-red fixed point near the scale of electroweak symmetry breaking, in the low TeV 
range, and natural values for the model parameters, the energy shifts in the low-lying muonic lead 
• transitions are typically of the order of a few times 0.1 eV to a few times 0.01 eV. The energy 

' level structure of the Unparticle Uehling shifts are inferred using general methods for the scalar and 

pseudo-scalar cases and it is shown that the two cases contribute to the energy shifts with the same 
' ^ ' sign. It is shown that this conclusion is not changed even when scale invariance is broken and is in 

' fact relatively insensitive to the scale at which it is broken. It is pointed out nevertheless that the 

estimated magnitude of the Unparticle Uehling shift (based on some natural values for the model 
parameters) is a factor of 1000 — 10000 below the discrepancy in QED/nuclear theory and precision 
muonic lead spectroscopy from about two decades ago. We briefly comment on scenarios where 
' (-H the Unparticle induced energy shift, if it exists, may be experimentally measurable. One possibility 

^ • in this direction is if the UV-sector, from which the Unparticle sector arises, has a large number 

I ' of fermions. Comments are also made on the possibility of further studying muonic atoms, as a 

Oh' probe for beyond-standard-model physics, in the context of forthcoming experiments, such as those 

probing lepton flavor violation through coherent muon-electron conversions. For completeness we 
explore some of the astrophysical and cosmological consequences of a fermiophobic scalar/pseudo- 
scalar Unparticle sector. In the fermiophobic context we also estimate a minimum value for the 
, conformal invariance breaking scale. 
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I. INTRODUCTION 

Recently there has been much interest in the possibiUty of a scale invariant hidden sector that couples to the standard 
model (SM) The operators in such a theory have been referred to as "Unparticles" to emphasize their generally 

non-integral scaling dimensions. In addition to an n-body phase space resembling that of a non-integral number of 
massless particles, the correlation functions of the Unparticle operators have interesting non-trivial phases The 
above properties make this sector distinct from other beyond-SM extensions. 

A model with the above properties was proposed by Banks and Zaks (BZ) 01 . Their theory was a vector-like SU (3) 
gauge theory with A^f fermions in the fundamental representation. It was found that for a particular range of Nf the (3 
function at lowest order is negative (/3i ( A^f ) < 0) while the contribution to the /3 function at the next order is positive 
{P2{Nf) > 0) in 

This would mean that as the theory flows to lower energies the small coupling constant g would grow until it hits an 
infra-red fixed point where 

Pig*) ^ 

The theory is scale-invariant below this scale and the description in terms of the Banks-Zaks fields at high energies is 
replaced by one in terms of composite particles of a strongly-coupled scale-invariant theory. These composite particles 
may be identified with the Unparticles [1]. We are interested in exploring the case of a scale invariant fermiophobic 
Unparticle sector that couples only with the massless SM gauge bosons, specifically with a substantial coupling only 
to the photon. As we shall see a fermiophobic sector might be able to avoid certain constraints compared to a 
fermiophilic sector. Also, independent of considerations in our study there is great interest in a fermiophobic scalar 
sector in the context of electroweak symmetry breaking and the Higgs mechanism ^^l. Our main focus will be on 
a fermiophobic sector that couples predominantly to photons and the effects of the induced oblique corrections on 
atomic energy levels in muonic atoms. 

If one assumes that the scale invariant Unparticle sector is also conformally invariant, then the scaling dimensions (j) 
of the gauge invariant primary Unparticle operators are tightly constrained by requirements of conformal invariance. 
For an operator in the (li, I2) Lorentz spin representation the constraints are [^-ISl 

j > ?i + ?2 + 2 - Si^xhfi 

These bounds translate to 

ju > 1 ; jui > 3/2 ; juv > 3 

for the scalar, fermion and vector Unparticle operators respectively. These conditions are referred to as Mack's 
unitarity criteria [3| • We will impose these constraints on the Unparticle operators we work with in the present study. 
Although it is technically possible for a quantum field theory to be scale invariant but not conformally-invariant, 
examples are rare [?[. 

The Unparticle operators in the limit of exact scale invariance may be considered as a sum over resonances having a 
continuous mass distribution with no mass gap [8| . The Kallen-Lehmann spectral density of the Unparticle operator 
takes the form [l[ 

p{s^) oc (sy-^ 

This implies that for j > 2 the theory becomes very UV sensitive and may cause singular behavior jQ] . This suggests 
that we only consider operators with j < 2. This is in conflict with the requirements of Mack's unitarity for primary, 
vector Unparticle operators which require j > 3 [!| . Usually in any Unparticle model with couplings to SM fields there 
must also be additional contact terms between SM operators as first pointed out in Q. They are generated when 
we integrate out the ultra- heavy UV fields (BZ fields for example). Without some fine-tuning these terms generally 
dominate over the SM- Unparticle couplings and render Unparticle phenomenology to be of lesser interest. In the case 
of scalar/pseudo-scalar Unparticles where the Mack's unitarity constraint is j > 1 and the scaling dimension may be 
very close to unity these contact interactions may nevertheless be of lesser importance 6] . 

Apart from simplicity the arguments above will be our motivation for considering scalar/pseudo-scalar Unparticles 
(U) in our study. In the subsequent analysis we will use U to label a generic scalar/pseudo-scalar Unparticle, O for 
an Unparticle scalar and O for an Unparticle pseudo-scalar. Also, we label an ordinary scalar by 0, pseudo-scalar by 
(f> and a generic scalar/pseudo-scalar operator by The energy scale of the infra-red fixed point (where the theory 
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becomes scale invariant) is usually denoted by Au and that at which scale invariance broken by the parameter /z. We 
will assume that the Unparticle scale is in the low TeV range near the scale of electroweak-symmetry breaking, i.e., 
^ V. Specifically we will interpret this to mean Ay £ 246 GeV, 0(1) TeV]. There are stringent constraints on 

pp ~> hf ^ 77 + X 

where ft,/ is a fermiophobic Higgs, from the CDF and DO collaborations at the Tevatron [l^. In our case the coupling 
of the fermiophobic scalar/pseudo-scalar sector to the fermions and the heavy gauge bosons W^) are assumed 
to be very small or close to zero. Also, the coupling to the SU{3)c colored massless gluons is assumed to be much 
smaller than the coupling to the photon. In this case the interaction scale may be in the Au ^ v range we consider 
above, in the low TeV regime, and still be consistent with the collider constraints [Tl'|. We will estimate a lower bound 
for the scale of /i in the fermiophobic Unparticle case later in section III. 

Based on all the above arguments we will therefore consider Unparticle scalars and pseudo-scalars coupling to 
photons with an interaction scale A^ ^ v and satisfying 

1 < < 2 (1) 

Now if one considers two representative fermiophobic and fermiophilic effective interaction terms of the form 

A-^UF^pF^P ■ A]f'U^^-i^-^ 

it may be argued based on scaling arguments and dimensional analysis (12j l that in general 

> A,.^ 

This seems to imply, for the relevant values of j, that the coupling to the gauge bosons may lead to larger effects 
than the corresponding fermiophilic case jl^ . Although the true niche for probing such effects may be in high energy 
colliders (see for example [1, [Tl| - [T7| and references therein) , especially when the interaction scale is much higher than 
0(1) TeV, it is still interesting to explore the effects of a fermiophobic coupling in low-energy experiments. In the 
fermiophilic case for instance one can put interesting bounds on the Unparticle sector from atomic parity violation [is} . 

We are specifically interested in fermiophobic Unparticle contributions to muonic atom transitions. Although in 
this case, unlike atomic parity violation in the fermiophilic case [3], no symmetry is violated there could be atomic 
systems or regions in parameter space where the effect may be measurable. Also, if the fermiophobic Unparticle sector 
is such that it has substantial couplings only to photons and no other gauge bosons, the muonic atom transitions can 
play a unique role in constraining it. 

We chose to study fermiophobic Unparticle contributions to muonic lead transitions. It will be seen in this specific 
case that for natural values of the model parameters the estimated Unparticle induced energy shifts in low-lying 
muonic-lead atomic transitions may be of the order of a few 0.1-0.01 which is comparable to the bound state 
QED corrections from light-by-light scattering and the fourth order Lamb shift to higher orbital angular momentum 
transitions respectively. We will comment on cases where this value may be enhanced or measurable, for instance 
when there is a very large fermion multiplicity in the UV sector. 

Another point is that since oblique corrections may potentially involve new "heavy" particles running in the loops 
we may expect them to be a probe for Unparticles within a wide range of n values. We will discuss this aspect, in the 
context of the Unparticle Uehling energy shifts, in more detail in section III. There is also the motivation that many of 
the constraints on the couplings of an Unparticle to SM fields [l^ are more relaxed in the case of a fermiophobic sector. 
We will discuss the case of broken scale invariance and the requirements for evading astrophysical and cosmological 
constraints later in our study. 

The scalar/pseudo-scalar Unparticle propagator in the limit of exact scale invariance is [l[ 

/ d4xe^«-(0|r(O^(x)O^(0))|0) = -^4^h(g2 + ze)]^-' (2) 
J 2 sm(j7r) 

where 

167r5/2r(j- -H/2) 
(27r)2jT(j - l)r(2j) 

with j being the scaling dimension of the Unparticle operator. Note that since j is in general a non-integer there 
is a non-trivial phase e*^^"-'^'^ associated with the Unparticle propagator. Also note that in the limit of j — )■ 1+ 
there is no singular behavior and we recover the ordinary scalar/pseudo-scalar propagator. The propagator coefficient 
ylj/(2 sin jtt) has unit magnitude at j = 1 and diverges very close to j = 2 where the theory is very UV sensitive. In 
the rest of the region its magnitude is below 1. 
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The case of broken scale invariance may be parametrized by introducing an effective mass gap /i in the spectral 
decomposition leading to the scalar/pseudo-scalar Unparticle propagator fwL flGj 

d'x e'^^0\T{Ou{x)Oum |0) = [- {q^ - /i' + (3) 

2 sin(j7r) 

fi may be thought of as the scale at which conformal invariance is effectively broken. 

Let the fermiophobic Unparticle fields be gauge singlets under the SM gauge group. Then the coupling of the scalar 
Unparticle to two photons may be incorporated [H by a term in the effective low-energy action 

St - I d'x-^OF,,F^'^ + ... (4) 
J 4 A:^ 

Here . . . denotes terms that are suppressed by higher powers of the relevant energy scale and which have been ignored. 

is a scale relevant to the 7 — — 7 coupling derived from the fundamental Unparticle scale Ky. The coefficient c 
is assumed to be an C(l) constant factor. Note that when the Unparticle sector is being generated from a UV theory 
(such as the BZ theory) the coefficient of the above operator goes like T] 



0{l) 



A;, 



Muv 



and is in general not of 0(1). Here Muv is the scale of the UV physics and o?uv is the dimension of the UV operator 
coupling to F^i/F^''. But if one assumes, for example, that the scalar Unparticle operators Oi are being generated 
by the confinement ( {^i^i) — > Oi) of fcrmions (^'i) from the UV sector (in this case duv — 3), conservatively, the 
"effective" coefficient 



i=l 



VMu 



could naturally be of 0(1). The reason is that theoretically the number of fermions Nf in the UV theory is permitted 
to be large and is constrained only by the requirement of conformal invariance at g* . For example 

^ > AT, > ^ 
2 - ^ 38 

in the BZ theory Q and in a recent Technicolor inspired SU{Nt) x SU{Nu) model by F. Sannino and R. Zwicky [23| 



7* + 2 ^ ■> ^ J* + 2' 

where the critical anomalous dimension satisfies the unitarity bound 7* < 2. Thus we will interpret the effective 
interaction in Eq. Q as modeling the effects of various possible Unparticle scalar operators (O « ^i)^ i^^ ^ 

semi-realistic model, and take the coefficent c ^ 0(1) without loss of generality. But it is nevertheless important to 
keep in mind that apart from notions of naturalness nothing excludes a larger value for the coupling, for instance if 
there is very large fermion multiplicity (ie. very large N^) in the hidden sector. QCD-like models with a possibly larg e 
number of colors {Nc) and fermion flavors (iVp) is not uncommon, for example, in many string-inspired models [21| . 
Similarly the coupling of a pseudo-scalar Unparticle to two photons may be modeled by 

S;f=[d'x-t-dF^^F,, + ... (5) 
J 4 A4 



where 



P —-fa paP 



Ay is a scale relevant to the 7 — O — 7 coupling and the constant b is assumed to be of 0(1). All the assumptions in the 
scalar case again hold here. The above may be compared to the chiral anomaly induced j—tt'^ — j Wess-Zumino-Wittcn 
coupling 

S.,, - / ^° F"" (6) 
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FIG. 1. Photon vacuum polarization by an Unparticle scalar or pseudo-scalar with an arbitrary scaling dimension j. The 
vertices are shown with blobs to indicate that they are effective couplings coming from an effective low-energy action of the 
form (gl) or (O. 



leading to 

valid in the soft pion limit (Here Nc is the number of 'colors', A;'' and Z"' are the four- momenta of the two photons) 
and the 7* — 7 — vertex (here is a pseudo-scalar meson like tt'^, ry or 77 ) 

T/x > -le I I e^upaP e'^q 

proposed in for = —q^ — >■ ±00 and recently considered in [2^ in the context of meson-photon transition form 
factors in the charmonium energy range. In the above expression g'^ is the four-momentum of the off-shell photon 7*, 
p'^ is the four-momentum of the pseudo-scalar meson ((>, is the polarization vector of the outgoing on-shell photon 
7 and /tt ~ 93 MeV is the pion constant. 

The corrections to muonic atom transitions due to photon vacuum polarization are induced by diagrams like those 
in Fig.[T] It is our aim to estimate the magnitude of such contributions to low-lying muonic atom levels. We mention, 
as an aside, that processes such as those in Fig. [T] for the case of low-energy QCD (where U{k) is now again a pseudo- 
scalar meson like tt" or ry) contribute for instance to muon (g-2). Using Chiral perturbation theory based on Eq. ([5]) 
and (jj vector-meson dominance, in the relevant range y/s < 0.6 GeV, a diagram such as that in Fig.[T] for the 7r° is 
expected to contribute 



,(7r%, ^^ < 0.6 GeV) = (0.13 ± 0.01) x 10 
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A similar contribution from the fermiophobic Unparticle sector, for typical model parameter values, is expected to be 
very small and within experimental limits, since the QED Kernel K{s) [25] has a steep cut-off around 0(1) GeV. 

Let us now proceed to calculate the vacuum polarization functions and estimate the induced effective Uehling 
potentials for various cases. 



II. OBLIQUE CORRECTIONS AND THE UNPARTICLE UEHLING POTENTIAL 

We are primarily interested in possible oblique corrections to the photon propagator due to scalar/pseudo-scalar 
Unparticles as shown in Fig.[T] If they exist we expect such vacuum polarizations to modify the photon propagator 
over the usual SM corrections. These scalar /pseudo-scalar oblique corrections can show up potentially as very tiny 
energy shifts in muonic atom energy levels or in the anomalous magnetic moment of the muon. As we have mentioned 
previously our main focus will be on possible corrections to the atomic transitions in muonic atoms. We will first 
consider the case of perfect scale invariance /i — ?> 0. The case when /J. 7^ will be discussed in section III after we have 
determined a minimum value for fi in the fermiophobic case as dictated by astrophysical and cosmological constraints. 

Using the Feynman rule for Eq. Q (see for example [2^) the scalar Unparticle contribution to the photon polar- 
ization tensor is 

^^^^^ c^A, r d^k [q.{q + k)g^^-q<-{q + kY] [g • (g + fc)g^ - g„ (g + fc)'-] 
2A^^' sin(j7r) i (27r)4 [(^ + /j)2 + [_(/.2 „ ^2) + jg]2-j 
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This simplifies to the the expression 



with 



where 



c2Aje''(2-J)'rr(3 - i) 
'2A^^'sin(j7r)r(2- j) 



dx dy 5{x + y — 1) ^ 



dH q^{x-lf + (PI2) 

(2^)4 [;2_ A(g2^^2)+j£]3-J 



A(g^/i2) = x{x~l)q^+ yfi^ 

We note that Eq. ([7]) has the expected gauge invariant structure and satisfies the Ward identities. 
For the pseudo-scalar Unparticle case the polarization tensor is 



' O - ~ 2A2J' sin(i^) J (2^)4 [[q + k) 



d*k [e^^'"'exc...qpiq + k)^q^{q + k)^] 



+ ie] [-{k^ - n'')+ie] 



2-i 



This may be simplified to give 



where 



■TT . 2 2 ■^ 6%e»(2-J>r(3-i) 
iU^iq^,^i^,j) = +--^J 



2A^^sin(j7r)r(2-j) 



dx dy 5{x + y ~ 1) y^ 



(/V2) 



(2^)4 



3-i 



and again 

A(g2,/x2) = x{x-l)q^ + yfi'^ 
Regularizing the momentum integral in Eq. ([5]) using dimensional continuation yields 



^o((7^/i^j) 



dx dy S{x + y — 1) y 



327r2A^^' sin(j7r) 

For the pseudo-scalar case, again regularizing via dimensional continuation 

^a(9^ j) - f dy S{x + y-l) y'-^ [^i^ A^-^ (g^, ^2) 

327r2A7 sm(j7r) J [ J U - 1) 

We will address the j — >■ l"*" limit of the expressions in Eqs. PT|) and (IT^ later. 
The Unparticle polarization functions may be renormalized as 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



(13) 



so that as q — >■ the residue of the photon propagator tends to unity. 

Note that the non-trivial Unparticle propagator phase e*(2-i)'r (joes not make an appearance in Eqs. (|TT|) and (|12l) . 
It is found from explicit calculation that the phase is removed during the evaluation of the Euclidean loop integrals. 
Therefore the Unparticle vacuum polarization tensor does not have a complex phase from the propagator contribution 
and any imaginary part that Iln(q^ , j) picks up subsequently, if at all, should come from the kinematic region that 
q^ occupies. 

We may likewise calculate the polarization functions for ordinary scalars {(j)) and pseudo-scalars (0) with a two- 
photon coupling. The ordinary pseudo-scalar in this case may be compared to an axion-like particle (a) coupling to 
photons 



I 
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where the scale fa and the pseudo-scalar mass nia are independent of each other. Substituting j = 1 in Eqs. ([8]) and 
(jlOp and integrating over the loop momenta using dimensional regularization we get after an MS subtraction 



^I(g^m|)= - 



dx dy S{x + y — 1) 
dx dy 5{x + y — 1) 



A(g^TO^) log —-^+q\x-iy log 



A ((7^m|) log- 



Af2 

A"(<72,m|) 



M2 



(14) 



Here 



A {q^,m1) = a;(x- 1) + ym^ 



and we have used the fact, from Eq. ([2]), that 



A, 



2 sin(j7r) 



^ - 1 



^ -1 



M is an arbitrary subtraction scale. It must be mentioned that this dependence on M will be removed when we 
explicitly introduce other higher order interactions in the effective Lagrangian. A natural choice is to take M ~ 
(A^, A^) ~ V, the energy scale of the relevant interaction. The choice would also ensure that in the m$ — ^ limit we 
can retain to good approximation, as far as numerical computations are concerned, just the lowest order terms in the 
effective Lagrangian. This is because with this choice, for the range of q we are interested in (specifically q < m^-), 
we have q^ \og{q'^ /PP) ^ q'^. For instance in the analogous case in QCD chiral perturbation theory (ChPT) a choice 
of M ~ Aqcd ^ 1 GeV is appropriate and one can keep to lowest order only the so called chiral logarithm term to 
estimate some of the low-energy QCD effects (see for example [13] and references cited therein). With this tacitly 
assumed we proceed to analyze the functional forms of the Uehling potentials in the cases of interest. 

If one is interested in bound states the modified electromagnetic four-potential A^, due to the vacuum polarizations, 
is given by 



AJ^) = 



q 



(27r)4 
d^q 



1 



il{q^) 

n(g2) 



(15) 



where i Gn^°*°" (q'^) is the photon propagator, J^ZuicSq) the source four-current in momentum space and A'"^\q) is 



the four- vector potential without any vacuum polarization corrections. As is well known 1 - 
constant for vacuum 



'il{q^) acts like a dielectric 



n(g2) 



£0 



Thus, the imaginary part of the polarization tensor corresponds to the vacuum becoming absorptive. In general 

fl(g2) = ns,,(g2) + Ui^{q^^^^J) + no„,„(g2) 

In our analysis we are interested in the corrections solely due to the Unparticle contribution Ilii{q'^, We noted 

previously that the non-trivial Unparticle propagator phase e^^^-i)'^ jg effectively cancelled during the evaluation of 
the loop integrals. From Eqs. pT|) . and (ITil) we observe that for 



x (x ■ 
x{x - 



l)g2 



9m [nuiq^^l^J)] ^ 








and there is a branch-cut starting at q^ — ji^ and q^ — m\ in analogy with — e~ vacuum polarization in QED. 
We will approximate an atomic nucleus of charge Z initially as a static, point-source. For a static source of the 
electromagnetic field the momentum transfer is space-like (q^ < 0) and the polarization tensor Iln{—q'^,^'^, j) is real 
which would imply that there is no absorption in vacuum under these circumstances. 
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U 



FIG. 2. Direct scalar/pseudo-scalar Unparticle exchange in a fermiophilic scenario where the Unparticles have substantial 
fermion couphngs. These possibihties are heavily constrained, especially when ^ — ^ 0, by searches for a 5"^ force, atomic parity 
violation and other experiments [T^. [T^. [2^ . 



Approximating an atomic nucleus of charge Z as a static, point-source the electromagnetic 4-current may be 
calculated in momentum space to be 

In our convention the charge e is intrinsically negative. Then the expression for the modified potential becomes 



Ze 



(27r)3 



(16) 



The I{u{<f 1 li) contributes to a Uehling potential (in a way similar to e+ — e vacuum polarization in QED) and 
gives 



Ze 



(27r)3 " 1 'r- / ='00 



(17) 



This is the correction to the electromagnetic potential due to fermiophobic, scalar/pseudo-scalar Unparticle vacuum 
polarizations. 

By simple dimensional analysis of our expressions in Eqs. (fTTj) and ([T2|) . using Eq. (|17p . we may expect the Unparticle 
corrections to the electromagnetic potential in the limit — > to be of a functional form 



1 



r2j + l 



(18) 



for j ^ 1. From Eq. (1141) for the ordinary scalar/pseudo-scalar case, in the limit m$ -> 0, we infer similarly that the 
potential may go dimensionally as This seems to suggest that, due to the generally non-integral values of the 

scaling dimension j, the fermiophobic Unparticle induced potential may have rich functional dependences. We will 
have more to say on this issue when we discuss the j — > 1 limit for Unparticle operators in the context of atomic 
energy-shifts. 

We can compare this fermiophobic scenario (where the correction to the potential is through photon vacuum 
polarizations) to the case of direct Unparticle scalar exchange between charged/uncharged fermions ( such a coupling 
could be generated radiatively at a higher order from fermiophobic couplings in the case of the charged particles) as 
shown in Fig.[2j In the direct exchange case (fermiophilic) 



(2^)4 



Therefore in this scenario we expect the correction to the potential to go like 



(19) 
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Of course in the case of broken scale invariance the direct exchange will be Yukawa suppressed (by the scale breaking 
parameter fj,) leading to a finite range for the force. We are merely interested in the behavior of the direct-exchange 
potential as a comparison to the vacuum polarization case. As discussed in section I we will require 1 < j < 2 to 
be consistent with the Mack's unitarity condition and the requirement of UV insensitivity. This immediately implies 
that for the fermiophobic case under the point nucleus approximation the potential is always more singular than 
This has some implications. We will look more closely at the sign of the Unparticle Uehling corrections later, but for 
now let us assume that as r — the potential energy is negative. Consider an e-ball, near r = 0, of radius r^. Then 
from the uncertainty principle the mean energy at is 



2j + l 



where is a positive constant that is determined from some explicit computation. We have (2j + 1) > 2Vj > 1 and 
this implies that for arbitrarily small r^, the mean energy (E) is arbitrarily negative for a negative potential energy. 
Therefore the particle must fall to r = signifying an instability (see for example [1^). This may also be reasoned 
semi-classically from the criterion that the negative potential energy must be less singular than the centrifugal term, 
l{l + l)/r'^ in the Hamiltonian, to have a stable bound state. Note that for a positive potential energy (^^ < 0) 
this constraint does not apply. Also note that if we had not assumed the scale invariant Unparticle sector to be also 
conformally invariant then j > 1 is no longer necessary, nevertheless for a negative potential energy the scenario is 
still problematic for Vj > 1/2. In this case even for j = 1/2 we have to require the condition 



e < 



8m 



to get a stable bound state as is familiar from quantum mechanics [29[. 

Thus it seems, at least preliminarily, that the correction to the electromagnetic potential in the direct exchange 
case is in general less singular than the fermiophobic case for any given value of the scaling dimension (j) under the 
point nucleus approximation. These expectations are modified to a large extent by some physical considerations. The 
first obvious modification comes from the fact that near r = the finite size of the atomic nucleus becomes important 
and we would have to calculate the correction to the electromagnetic potential with an appropriate nuclear density 
profile. This would soften the singular nature of the potential near r — 0. One may consider the expression in Eq. (jl8p 
as being roughly valid far away from the atomic nucleus in which case it may be approximated by a point source. 

More importantly one needs to be careful in the evaluation to get consistent energy shifts as j — 1+. From Eqs. pT|) 
and ([T^ we see that the j 1+ limit is not continuous. It seems to imply that if one calculates a potential and 
an energy shift using the quoted Unparticle polarization functions the energy shift increases without bound as one 
approaches j = 1. This does not seem physically sensible. In the j — 1+ limit, maintaining continuity throughout, 
one must expect to recover the ordinary scalar (0) and pseudo-scalar ((/>) situations. This feature seems to be peculiar 
to the Unparticle vacuum polarization diagram we are considering and is related to the fact that the momentum 
integrals diverge as A^-' , where A is some momentum cut-off, and are only quadratically divergent when j = 1. In a 
similar calculation for the fermiophilic vector Unparticle contributing to muon (g — 2) (where the photon in the loop 
is now replaced with a vector Unparticle) the relevant momentum integral, for /i ^ 0, is [isj 



m. 



dH 



1 



P - Av(g2, ml_)+ie 



4-j 



A-'" 



m,, 



-) 



m^- is the muon mass and A^ is the interaction scale for the fermiophilic coupling. The limit j — )■ 1+ is now 



continuous and one obtains 



mf,_ A- 



167r2Av((z2, 



as expected for an ordinary vector particle in the loop (like the photon). 

So to make the expressions (jll[) and (1121) consistent with the j — > 1^ limit and maintain continuity we will consider 



A^-l(cz^M') 



J 



1 



(M^y-' loe 



A(g^A^^ 
M2 



(20) 



where a term that would lead to an infinite energy shift in the j — > 1"'" limit has been subtracted out. Here M is again 
an arbitrary scale. As we mentioned previously for the ordinary scalar/pseudo-scalar case a natural and appropriate 
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choice is to take M ^ Ay, the scale of the infra-red fixed point, where the theory becomes scale- invariant and the 
description is in terms of Unparticles. The crude prescription in Eq. (|20l) is the one that we will follow in the present 
study to estimate Unparticle Uehling shifts for values of j away from 1 and make the correspondence with the j — 1 
case. Note that for j = 1 the prescription gives the ordinary scalar/pseudo-scalar expression. We will be interested 
mostly in the case when j is close to 1 where, as we shall see in the next section, the Unparticle induced atomic energy 
shifts have their maximum values. Also, as mentioned in section I, very close to j = 1 the SM-SM contact terms are 
relatively less important numerically compared to the case when j is very far from unity [6]. We may therefore be 
optimistic that the above prescription will certainly capture the essential features of the Unparticle Uehling shifts in 
the interesting region near j = 1. 

In the /i — >■ limit we may factor out the q dependent terms in Eqs. pT|) and (IT^ . Performing the integration over 
the Feynman variables we may show that the polarization functions for the scalar and pseudo-scalar Unparticles are 
exactly equal y j: 

^o(-g^o,j) = n<5(-q2,o,j) 

assuming b = c and A-y = A-y. For the ordinary scalar/pseudo-scalar cases this equality may be understood readily 
in terms of the optical theorem. One consequence of our approximation in Eq. ()20|) is that this relation is no longer 
exact when j is sufficiently far from unity. We will therefore retain both the scalar and pseudo-scalar expressions in 
all our analyses. 

Let us now perform the analysis incorporating the finite extent of the nucleus. Assuming a static situation we have 
Let 

jL..M=P{r}= -Zef{r) 

where /(r) is a suitable function, which we assume to be spherically symmetric for simplicity, describing the nuclear 
charge density profile. It must be suitably normalized such that 

j d^rf{r)^ 1 

With this choice one may now perform the integration over the angular variables to obtain 

Mo(f) = Vu{r) -ZeJ ^ ^^^jMll II^(-g^^^,) fi\q]) (21) 

where 

= J dVe'^- /M (22) 
We will choose a simple gaussian nuclear charge density profile 

/M - /o e-'^'/(20 (23) 

with 



(27rC)3/2 

to give the correct normalization. Note that at r = the nuclear charge density is 1/e times its value at the 
origin. The gaussian nuclear profile may be considered to be an approximation to a more realistic two-parameter 
Fermi distribution. Usually the gaussian density profile is considered more appropriate for light nuclei ^30,] . In the 
next section we will see that we are interested in studying muonic lead, which is not a light nucleus. Nevertheless 
even with the simplistic choice of a gaussian profile we may expect the analysis to capture all the main features of 
the finite lead nucleus. With this ansatz, for the nuclear profile, we expect any deviation from the actual case to be 
of 0(1) and hence comparable to our ignorance about the factors c, b in the effective Lagrangian. The main reason 
for choosing a gaussian profile is to simplify the subsequent analysis. The choice would also help us glean simpler 
analytic results which otherwise would have been more intractable and subject to numerical methods solely. 
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As previously mentioned for now we will assume exact scale invariance and set ^ — >■ 0. We will comment on this 
aspect later in section III and will incorporate a non-zero /i at that time. It will be seen that the main conclusions 
are not drastically changed for the case of broken scale invariance. Defining the dimensionless variable z = |g| r we 
may write Eq. (|21l) for the two cases, using Eqs. (ITT|) and p^ . as 



Zec^ A, 



647r4A^-'' sin(j7r) jq 
(M2p-i 



dxdyS{x + y- 1)?/^ ^(1 - a;) [{1 + j)x - j] 



Zeb^Aj 
647r4A^^ sin(j7r) Jo 



dz z sin z e 
1 



log(x(l - x)) + log 



dx dy S{x + y — 1) [x{l — x)y 



(24) 



(25) 



dz z sinz e '■^ Z'^'" ' 



log(a;(l - x)) + log 



One may now perform the integration over the dimensionless parameter z and the Feynman variables. The result may 
be expressed in terms of the exponential function and the generalized hypergeometric function F^[{ai, .., Cp), (6i, .., hq); w] 
in a compact form 



Voir-, 3)= - 



Zee-" 


BiU) 


r2 

e 2C 


^ S2(j) 


^2 


64 7r4 


(3/2 ^ 


(5/2 




BiU) 






j,2 g-2c 


64 7r4 


(3/2 ^ 


(5/2 



[(l,l),(2,5/2);rV(2C)] 
F^ [(l,l),(2,5/2);rV(2C)] 



(26) 



Bi{j), B2{j), Biij) and B2{j) are functions of the scaling dimension j and are given by 



V^A, (M ^) 
V2A^^' 



2V'-1 



sm jTT 



2V'-1 



2 + j(i - 2)(2j - 7) - (j - 4)(j - 3)(1 + j(j - 3)) H[3 - j] 



B2U) 
Biij) 

B2U) 



3 V2 A^"* sin j'tt 
\f2Ki^ sin jTT 

^A,[Apy-' 

3 V2 A£"' sin ITT 



- 3j + 1 
.7 (j - 4) {3 - 3) 

-.?' + 3j-l ' 



3 U - 4)2 (j - 3)2 
log (2M2(e^^-2 



(27) 



.J0--4) (j-3)_ 

-5+j(5-j) + (j-4)(j-3)H[4-j] 
J (j - 4)2 (j - 3)2 



1 



j(.?-4)(.7-3) 



log (2M2(e''^" 



(28) 



J (J--4) (j-3) 



where is the Euler-Mascheroni constant and II[n] is the n"' harmonic number given by X]fc=i(l/^) ^^r integral 
values of n and by the Eulcr integral representation 



dx 





for n ^ +Z. Note that -\/( has dimensions of length and hence the above expressions for the potentials have the 
expected dimensions. The functional forms of the potential energies in the pseudo-scalar Unparticle case, eV^{r,j), 
for two different values of the scaling dimension j are shown in Fig.|3l The profile indicates that for the pseudo-scalar 
Unparticle the potential energy is negative. It is seen that the singular nature of the Unparticle Uehling correction, 
near r — 0, has been eliminated by taking into account the finite size of the nucleus. We will explore in more detail 
the sign of the Unparticle potential energies shortly. 

The Unparticle Uehling potential would cause a small energy shift in the atomic energy levels. The energy level 
structure in a bound system may be deduced (without explicitly calculating all the energies) by studying the relevant 
potential We will apply some of these methods to study the level structure of the energy- shifts (5E^) due to 
Unparticle Uehling potentials and ([25)) . 
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FIG. 3. Potential energies in muonic lead for the pseudo-scalar Unparticle case based on Eq. (|25p . assuming perfect scale 



invariance ^ 0, for j — 1.01 (dashed line) and 1.15 (small-dashed line). We have taken 6 ~ 0(1), — 4fm (to model the 
case of a pT ~ Pbg^** nucleus) and adopted a reference value for the interaction scale A.y ~ 246 GeV. Note that the potential 
goes to zero very rapidly away from r = and has a constant value at the origin implying that there is no singular behavior. 



A general result 32] is that for a potential V(r) depending on 
the energy levels are ordered as 

E(n, I) I E{n, I + I) 

Here n is the principal quantum number and I is the angular momentum quantum number. Note that for the simple 
Coulomb potential Vii^) — for non-zero r and gives the familiar result that energy levels depend only on the 
principal quantum number n. For the Unparticle Uehling potentials being considered we see that the situation may 
be complicated somewhat by the fact that for the finite nucleus the Laplacian of the potential energy may change 
sign. Thus strictly speaking we must consider if 

(eV^Vu{r,j, i?)) I (29) 

Here the averaging (. . ■)n,i+i is over the relevant (n, / -|- 1) wavefunction [s^]- Depending on the sign of the above 
expression we would have a relation between the Unparticle induced energy shifts 

Let us proceed to study the Laplacian of the potentials in Eqs. and (|25l) . As mentioned before we will consider 
the case of /i ~ presently and concentrate on the low-lying states, i.e. states near r = 0. It is noted, from Eqs. (P7)) 
and (|28p. that for both the scalar and pseudo-scalar Unparticle potentials 

S2(j) > 
Bi(j), Ml) > 



for Vj G [1,2). It is also seen that 



BiU) » B2{j) 
BiU) » Mj) 
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Thus we infer that both pseudo-scalar and scalar Unparticle Uehling potential energies are negative Vj € [1,2). This 
reaffirms the indications from Fig.|31 It also means, as indicated by the i?i(j) and i?i(j) coefficients, that the first 
term in the Unparticle Uehling potential dominates over the second term. The Laplacians of the Uehling potentials 
in Eq. (|26p are given by 



eV^Vo{r,j,0) 



Ze^ 
64 



B,{j) ^ (r2 - 3C) + B2{j)^\ 12re ~ 3e€ V2^C'/' Erf 



(V2 



^/K 



(30) 



+ r3 (r2-3C) n [{l,l},(2,5/2);rV(2C)] 



eV^V^{r,j,Q) =- 



64 7r4 



Bi(j)^ {r' ~3C) + B,U)^[l2re - 3e^V2^ef'Eri 
+ r' (r2-3C) Fl [{1, 1}, (2, 5/2); rV(2C)] 



e 2c 



where Erf [z] is the error function defined as 



2 

Erf[z] ^ — dt 
V"' Jo 



The Laplacians of the Unparticle Uehling potentials correspond to vacuum charge densities created by the Unparticles 

Po{rJ)= - eoWo(r,j,0) 
Pl^ir,j)= - eoWg(r,j-,0) 

Since the potentials in Eq. fall off faster than 1/r (as r — > oo) the net, Unparticle induced, vacuum charge vanishes 
trivially 



Qir - / d'rp^-{r,j) = I d\ [-eoW;,(r,j, 0)] = -eo Ymi^j^da-VVu{r, 



i,0)- 



as we should expect. 

Now from our analytical expressions it is observed, as mentioned above, that for typical values of the scaling 
dimension j and radial coordinate r, the first term in the Unparticle Uehling potential dominates. It is seen from 
Eq. (|30p that the first term changes sign at r = \/iC,- Thus we may suspect that the sign of the Laplacian for the 
complete Uehling potential would also be dominated by the sign of the first term. Therefore we have for both the 
potentials 



eV2Vb(r,j,0) ^0 V r ^ 
eV^V^{r,3,Q) ^0 V r ^ \/3C 



(31) 



These expectations are confirmed by explicit numerical computations where it is found that the Laplacian of the 
Unparticle Uehling potential indeed changes sign at \/iC,. 

We are interested in the low-lying atomic states which will have relatively large energy shifts due to the Unparticle 
Uehling corrections. The higher angular momentum {p, d, . . .) wavefunctions are spread out from the origin. Specif- 
ically, it is seen that they are mainly non-zero for r > %/3C for semi-realistic values of C that will model an atomic 
nucleus ( for example C ~ 4fm^ for the lead nucleus). From Eq. ((29| we therefore note that due to this, the sign of 
( Vu{r^ j, /^^))n, i-i-i will be dominated by the sign of the Laplacian in the region r > \/3C. Thus for the Unparticle 
scalar and pseudo-scalar Uehling potentials 



< 

< 



(eV^Vo{r,J,^i^)) 

\ / n, (+1 

(eW'Vg{r,j, fi^)) 

\ ^ / n,l+l 

This implies that in the case of perfect scale invariance we have 

{n, I, j, 0) < SE'^{n, I + 1, j, 0) 
SE^{n, I, j, 0) < SE^{n, I + 1, j, 0) 



(32) 



(33) 
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This would mean that if we were looking specifically at ^ = 0, 1 and 2 states near r = the Unparticle Uehling 
corrections would be ordered as 

SE^^^U) < 5EZ,{j) < 5El,{j) 

We will explore the level structure of the Unparticle Uehling shifts in the next section and confirm explicitly the 
general results above. 



III. MUONIC ATOMS AND ENERGY SHIFTS 

Muonic atoms are created when muons are stopped in matter [s^ . The muon initially undergoes scattering and as 
it loses energy it may be captured by one of the atoms in a higher orbital. From here it cascades down, via various 
processes, to the innermost orbits where it persists until decay. The whole cascade process is expected to occur within 
~ 10~^ — 10~^^ s and the muon spends ^ 10^^ — 10^^ s in the inner orbits until decay. During cascade and its time 
in the innermost orbits we may probe the energy levels of the system (see (33 | and references therein for instance) . 

Muonic atoms are especially suited for probing the effects of oblique corrections to the photon propagator [s^, [sj] . 
In the QED induced hydrogenic Lamb shift [35] the 2si/2 state is found to lie above the 2p]^/2 state by 

<5£;^:r, = E2S,,, - E2p,,, ~ 1058 MHz 

which corresponds to about 4.4 x 10~^eV. The QED vacuum polarization contribution to the Lamb shift fsoj on the 
other hand is found to be 

6E^:^, = E2S,,, - E2p,,, ~ - 27.1 MHz 

For later comparison to the Unparticle case we note that the above — e~ vacuum polarization contribution is about 
1 X 10^^ eV. Thus in ordinary atoms the vacuum polarization effects (oblique corrections) contribute less significantly 
than the other QED corrections. Things are dramatically different in the case of muonic atoms owing to the much 
higher mass of the muon compared to the electron. Now, the energy shift due to vacuum polarization may be related 
to the Uehling correction dV{r) as 

<5i?;r= J d'rKiir)eSV{r)%Ar^ (34) 

We are specifically interested in the low-lying atomic states and the S and P states are particularly attractive since 
the modification to the potential, as we saw previously, mostly occurs near r = and falls off sharply as one goes 
away from the center. Taking the wavefunction 5'„/(r^, to lowest order, to be the Schrodinger wavefunction we note 
that near r = 

Thus we expect that in muonic atoms, where an electron in the inner orbital is replaced by a muon, the energy shift 
due to oblique corrections would be enhanced as 

Putting the masses in, (m^-/me-) ~ 10^. This is a huge enhancement to the vacuum polarization contribution. 
In fact it is found that in muonic hydrogen the 2si/2 state lies below the 2pi/2 level with QED vacuum polarization 
contributing about 206 meV compared to about 0.6 meV from other QED corrections '34'] . Another way to understand 
this enhancement is by looking at the wavefunction profiles in muonic atoms and noting that the wavefunction 
penetrates into the nucleus much more than in the electron case. In all our numerical computations we use the 
hydrogenic Schrodinger wavefunctions to lowest order, with m^- replaced everywhere with m^- and taking the 
appropriate value for the atomic number Z. One must in principle use the bound state solutions to the Dirac 
equation but for the purposes of our estimations the Schrodinger wavefunctions are sufficent. One aspect that we 
would like to point out though is that the bound state solutions to the Dirac equation are more localized near r — 
compared to the Schrodinger wavefunctions and hence will in general yield a slightly higher estimate for the corrections 
to the energy levels. Any error that comes from this simplification will be at most an 0(1) factor and similar to our 
ignorance regarding the coefficients in the Lagrangian or the interaction scales. Thus all calculated energy shifts in 
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our study must be interpreted as accurate only up to undetermined 0(1) factors to accommodate our ignorance of 
the interaction scales, coefficients and other approximations we have made. 
The typical binding energy of an atom with atomic number Z goes as 

But from Eq. we see that 

since 

Based on the above expressions we expect 



z' 



This motivates the reason for choosing a high Z system while probing for the existence of tiny Unparticle induced 
energy shifts. 

Based on the above considerations two very favorable systems for probing Unparticle vacuum polarization induced 
energy shifts are muonic mercury — HggQ°) and muonic lead — Pbg§^). We will choose as our reference system 
muonic lead. With Zp,-, = 82 and Ap^ = 208, some of the relevant scales in the muonic lead system are 

m^- ~ 0.1 GeV, A^:.„^,„„ ~ 12 fm, r^;,^ ~ 3fm , Cciou= « i?o 4f - 7fm 

where -Rq — 1-2 fm. From these crude estimates we see that the muonic Bohr radius in lead is much smaller than the 
extent of the lead nucleus or the muonic Compton wavelength and hence the muonic wavefunction will penetrate into 
the nucleus to a large extent compared to electronic wavefunctions. This results in an enhancement of the Uehling 
energy shift as we argued previously. It also implies that finite nuclear effects would be much more important in the 
low- lying states of muonic atoms, especially the muonic lead system we are considering js?!. l38l|. 

Another point to note is that for lead, with Z — 82, the parameter aZ used in perturbative QED calculations is 
no longer small. This makes the QED calculations more complicated and leads to larger theoretical uncertainties. 
Table I shows some of the QED corrections to the low- lying states of muonic lead. The interested reader is referred 
to theoretical details in fSJl and references therein. 



TABLE I. QED corrections, in eV, for the low- lying states in ii — Pbg2^ (see [13, and references therein). 



QED corrections 




2Sl/2 


2Pl/2 


2p3/2 


3Pl/2 


3^3/ 2 


3^3/2 


3^5 


Electronic Uehling and Kallen-Sabry 


-67864 


-19537 


-32648 


-30082 


-10871 


-10334 


-10605 


-9941 


Electronic Wichman-KroU 


492 


244 


348 


335 


160 


160 


186 


180 


Muonic Uehling corrections 


-248 


-43 


-45 


-34 


-14 


-11 


-1 


-1 


Leading self energy corrections 


3220 


696 


348 


649 


149 


224 


-44 


51 


Higher order self energy corrections 


153 


25 


65 


58 


21 


20 


8 


6 


Electron screening 


-5 


-25 


-13 


-13 


-52 


-54 


-37 


-39 


Recoil correction 


-382 


-87 


-111 


-95 


-30 


-26 


-15 


-14 



Similarly the various finite nuclear effects can be quantified to a large extent in muonic Phif (see M.-y. Chen 
(1970) [371], H. F. Skardhamar (1970) [381 and A. Haga et al. (2007) [H, for example). In Table II we give the current 
estimates of nuclear polarization (NP) effects in — Pbg2^ for three different models of transition densities as quoted 
in [11]. 

TABLE II. Estimates of the NP effects (eV) in /i" - -Pbg"** taken fr om [39|. The gauge dependences are shown in 
brackets. 
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Nuclear polarization corrections lsi/2 


2Sl/2 


2Pl/2 


2p3/2 


3Pl/2 


3P3/2 


3^3/2 


3^5/2 


TGT model -2727(4) 
RIN model -3599(10) 
JS model -5721(28) 


-463(1) 
-611(4) 
-930(8) 


-1357(7) 
-1590(10) 
-2178(13) 


-1425(9) 
-1656(10) 
-2214(7) 


-561(4) 
-690(3) 
-929(3) 


-749(1) 
-914(1) 
-1179(2) 


-226(0) 
-239(0) 
-280(0) 


-43(0) 
-42(0) 
-38(0) 



It was noted in the pas t that there is a discrepancy in the A2p and A3p NP calculations with results from muonic 
lead spectroscopy [iol l4l|. This discrepancy seems to have been partially tackled in the work of A. Haga et al. [s^ 
and the current discrepancy for A2p is in the ball-park of about 50 eV as shown in Table III. 

TABLE III. Comparison of Ap splittings (KeV) in /i" - Pb^§^ |39|] based on the QED results in Table I and NP 
calculations in Table II. 



^l- - Pb 


82* level 


TGT 


RIN JS Exp. 


2p3/2 - 

3p3/2 - 


- 2pi/2 
' 3pi/2 


184.858 
47.231 


184.846 184.829 184.788(27) 
47.208 47.225 47.197(45) 



Historically the spectroscopic information from low-lying muonic lead transitions was used to constrain parameters 
in the theoretical nuclear calculations. In Table IV the results of precision measurements on some of the low-lying 
transitions in muonic lead are shown. It is noted that the experimental uncertainties are typically of the order of a 
few tens of eV. 



TABLE IV. Precision measurements on some of the low-lying /i 



PboS^ transitions 14 



Muonic 
2p3/2 
2pi/2 



3d: 
3cij 
3d. 



3/2 
5/2 
3/2 



3P3/2 
3Pl/2 

2si/2 

2Sl/2 



transition 

^ lSl/2 
^ lSl/2 
^ 2pi/2 
^ 2p3/2 
^ 2p3/2 
O 2si/2 
O 2si/2 
O 2pi/2 
^ 2p3/2 



Energy'^ (KeV) 

5962.770(420) 

5777.910(400) 

2642.110(60) 

2500.330(60 ) 

2457.200(200) 

1507.480(260) 

1215.430(260) 
1030.440(170) 



° D. Kessler et al. (1975) 
P. Bergam et al. (1988) 



Energy" (KeV) 

5962.854(90) 

5778.058(100) 

2642.332(30) 

2500.590(60) 

2457.569(70) 

1507.754(70) 

1460.558(32) 

1215.330(30) 

1030.543(27) 



As suggested by Fig. [3] the states IS* and 25* are most affected by the Unparticle Uehling potential. A transition 
from 15' to 25* is nevertheless forbidden due to the electric dipole selection rule AZ = ±1. That leaves the possibility 
of a transition to an / = 1 state, that might still be sensitive to the Unparticle Uehling effect. We will consider, 
as the prototypical low-lying muonic transition, excitations between 15 and 2P. Based on a direct fit to muonic 
lead transition energy data (including higher level transitions) , keeping the variables in the Fermi distribution as free 
parameters, one may obtain an estimate for the discrepancy between theory and experiment in the 15 — 2P transition 



in -Pbgf [41[ as 



A^;^,,,. [lsi/2 - 2pi/2] - A£;c.,,./p., [lsi/2 - 2pi/2] ~ 227 eV (35) 
AE^^^_ [Lsi/2 - 2p3/2] - A£;c.:../p., [lsi/2 - 2p3/2] ~ -89 eV 



This estimate is a very conservative one based on a direct fit [4l| with a very poor x^/d.o.f = 187. We will take the 
above estimate as a crude measure of the discrepancy between measurement and calculation. It should be noted also 
that the NP values used in Eq. ([55)1 are obtained keeping the variables in the Fermi distribution as free parameters 
and are therefore not completely theoretical. 
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We will see shortly that, for typical values of the model parameters, this makes the identification of any possible 
Unparticle Uehling corrections in fi^ — Pbgj^ extremely difficult since the magnitude of any such correction will be 
found to be well below the NP uncertainties. This would imply that any Unparticle Uehling correction if present 
may not probably be unambiguously seen due to the data-fitting procedure required to extract information for the 
theoretical nuclear calculation in /i" - -Pb8^^ unless we can obtain such information independently. 

There are two scenarios where these conclusions could get modified. The first, as we commented earlier in section 
I, is if the UV sector has a very large fermion multiplicity (large Nf) which increases the total contribution from the 
fermiophobic Unparticle sector. The second is muonic atoms with intermediate- Z which may be more amenable to 
nuclear polarization calculations. 

We now proceed to analyze the IS" — 2P level corrections in various cases of interest. 

Let us first consider the case of ordinary scalars ((/)) and pseudo-scalars {4>). The relevant corrections to the 
electromagnetic potential is obtained from expressions in Eq. (ITi)) . In this case we find that the energy shift due to 
scalar oblique corrections is 



SEt, = 



where 



d'r |*t:;(r-)|2 eSV^r) 

dV|vi/^;(f)|2 



(36) 



167r2A2 



(2^ 



e''~-'{q')-'m / dx\A (q\ml) - q^x - I)' log 



A(_gVr4) 
M2 



A (g^ TO^) = x{x - 1) q2 + (1 _ x) ml 



and is the mass of the scalar as before. 
For the pseudo-scalar case similarly 



5Et 



4> 



_fq_ 



e'^-'~{fr'~f{q) 



dxA (g^TOp log 



A ((Z^mp 
M2 



(37) 



and 



A (q^ml) = x{x~l)q^+ (l^x)', 



denotes the mass of the pseudo-scalar. Let us consider the case when = and = 0. With this choice the 
above expressions may be evaluated numerically to calculate the energy shift in the muonic lead 15* — 2P transition. 
We take a, a ~ 2 and adopt the reference value A$ ^ 246 GeV with a renormalization scale M ~ A$. With these 
parameters the magnitude of the energy shift in the 15 — 2P transition for m$ = is estimated to be 



^\6Et,~SEt\^ \6Ei~6Ei\c, 0.140 eV^ O(0.1)eV 



(38) 



It is found that the correction to the transition, from the SM value, is positive for both scalars and pseudo-scalars. 
The above estimate may be enhanced due to a large multiplicity in the UV fermion sector and is not uncommon in 
many string-inspired QCD-like models [2l| as we commented earlier. 

Let us consider the same calculation for an Unparticle scalar/pseudo-scalar of scaling dimension j. We parametrize 
j as a deviation from the ordinary scalar/pseudo-scalar case by putting 



J 



1 



As argued in section I we require < 77 < 1. Then the energy shifts due to the scalar and pseudo-scalar Unparticle 
oblique corrections are respectively 



seZU) 



2A^^ sin(j7r) 
(47r)2 



d^q 



(27r)3 

-q-'{x-ir\ogA{q^f,')/M' + A{q',pi^) 



e'^-^\q^)-^~f{q) J dxdyd{x + y~l) y'^^ 
logA(g2,/i2)/M2" 



(39) 



(1 + 77) 



SE^iU) 



2Kij> sin(j7r) 



- nl 



d^q 



e'^-^q^)-^~f{q) / dx dy 5{x + y-l)y 



{M__ 

(47r)2 



A(9^M 



2 2, logA{q^^l^)/^'P 



(1 + 77) 



(40) 
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FIG. 4. Estimates on the magnitude of the energy shift IdE'ip — SEfs] due to possible scalar /pseudo-scalar Unparticie vacuum 
polarizations. We have taken jj, ~ 0, & = c ~ C(l) and A^., A^. ~ 246 GeV as the reference value at the renormalization scale 
M ~ An ~ A^, A^. As we have mentioned previously the quoted values for the Uehling shifts should be interpreted as being 
accurate only up to C(l) factors to accommodate suitable ranges for the interaction scales and other approximations. It is 
clear from the plot that the Uehling shifts from Eqs. H39p and (|4Up can under very general scenarios be in the ball-park of a few 
times 0.1 eV. 



where A(g^,//^) and A{q'^,fi'^) are as defined earlier in section II. 

The energy shift in the 15 — 2P transition in muonic lead may now be calculated as 

^Eg^is = SEg - SEE 

using Eqs. ([55]) and ([^0]) with the nuclear profile of Eq. ([25)1 . The results of the numerical computation, for /i — > 0, are 
shown in Fig.|4l We observe from the plot that the energy shift due to pseudo-scalar Unparticie vacuum polarization 
is generally smaller in magnitude than the corresponding scalar case assuming other parameters are identical. This 
is probably an artifact of our approximation in Eq. (|20p as we commented previously. As one approaches j — > 1^ 
though the energy shifts due to scalar and pseudo-scalar cases start to approach the same value of A£:2p_is - 0.14 eV 
as expected. Note also that in the j — ^ l''^ limit the ordinary scalar/pseudo-scalar case is recovered in a continuous 
manner. 

Therefore, it may be claimed that for /i — ^ the typical shift in the IS — 2P muonic lead transition due to the 
Unparticie Uehling potential is in the 0(0.1) eV range with reasonable assumptions about the parameters in the 
theory. This energy shift, if it exists, is similar in magnitude to QED corrections from the virtual Delbriick effect 
(light- by-light scattering), to higher angular momentum transitions, in /i^ — Pbg2* (3^ . For example the virtual 
Delbriick effect contributes 

AE^C.g^ 0.4 eV 

in /i^ — Pbg2*. The Unparticie induced energy shift may also be compared to the vacuum polarization contribution in 
the hydrogenic Lamb shift which is in the lO"*" eV range (ordinary hydrogen) or eV range (muonic hydrogen). From 
psp we see nevertheless that the discrepancy between measurement and theory of the 2pi/2 — lsi/2 and 2p^/2 ~ ^-^1/2 
muonic lead transitions is in the range of many eV j4Q . i41i] . This means that the estimated energy shift is about 
10^ — lO'^ times smaller than the discrepancies quoted in (|35|. We will analyze a more realistic scenario with /i 7^ 
later and will note then that the main features are not drastically modified from the simplistic /i — > case. 

It is also seen from the plot that the energy-shift is very nearly linear, with respect to the scaling dimension j, for 
values of j close to 1. In fact it may be verified analytically from Eqs. (|37P . (|39p and ([^01) that for j = l + rj with 
ry ^ 1 we may expand 

5Et, - SE^^iij) ^ Ci 77 + C2 + C3 + . . . 
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where dimensionful factors have been omitted and C„ are numerical coefficients. This imphes that very close to j = 1 
we have to leading order 

as suggested by the plot and Eq. (|20p . The above expression indicates that close to j = 1 the additional energy 
shift in the IS* — 2P transition compared to the ordinary scalar/pseudo-scalar case comes from the fractional scaling 
dimension part (ry) of the Unparticle. This seems very satisfactory in the sense that any correction to the ordinary 
scalar/pseudo-scalar case is just proportional to the "Unscalarness" of the Unparticle. This seems to be in the spirit 
of deconstruction [8] and is similar to the observation in [l^ 

'Tw(j) = i2-j)a4j^l) 

for the Unparticle production cross-section when the Unparticle sector is gauged under the SM. 

One also observes from Fig.|4] the onset of singular behavior as j — 2 as discussed in section 1. As pointed out this 
is a pathology arising from the fact that close to j = 2 the model becomes more and more UV sensitive. One may 
try to mitigate this singular behavior near j — 2 hy adding local contact terms Q ■ 

Concerning the sign of the Unparticle energy shifts it is inferred that there are no differences between the scalar and 
pseudo-scalar Unparticles. It found from our expressions in ([39)) and pO)) that the scalar as well as the pseudo-scalar 
Unparticle corrections to the 15 — 2P transitions are positive 

Ai?4-t ^ AE.l^U, + I A£; 1 (41) 

In the present case, for fi ^ 0, from Eqs. (|32|) and ([33]) we would expect the level structure of the Uehling shifts to 
be ordered as 

These expectations from Eqs. (1321) and p3p are confirmed by numerical calculations. The variations between the 
I = 0, 1, and 2 states are generally found to be of the order of a few 0.01 meV for n = 3 and of the order of a few 
meV for n = 2. 

Let us now turn our attention to the case of broken scale invariance where fJ- ^ 0- Even in the case of a fermiophobic 
Unparticle an effective interaction term with Higgs fields of the form 

is expected to be present. Here H stands for a generic Higgs field. As pointed out in the literature this interaction term 
is unique [isl l43j in the sense that it is a super-renormalizable term (since 1 < j < 2 from Mack's unitarity). Even if 
one starts with a near zero coupling of this term, at a high energy, it will be generated through renormalization group 
flow and will break the conformal invariance of the Unparticle sector, at a scale /i, when the Higgs develops a vacuum 
expectation value. Thus there is a very strong theoretical reason to expect scale invariance, in the fermiophobic 
Unparticle sector, to be broken at some scale after electroweak symmetry breaking. 

Moreover, from the observational viewpoint, it is expected that there should be very stringent constraints on 
massless/light scalar degrees of freedom from cosmology and astrophysics. For a scalar/pseudo-scalar Unparticle 
coupling to two photons with u, — )■ (i.e. perfect scale invariance) there are very strong bounds on the couplings from 
supernovae cooling [l9l |45| for example. Also, in the case of perfect scale invariance, there are very tight bounds 
on the couplings from big-bang nucleosynthesis (BBN) (45l |46|. These constraints effectively render any collider or 
low-energy experiments ineffective in probing the Unparticle sector. 

We would therefore like to estimate a lower bound on the scale invariance breaking effective mass (/i), for the 
fermiophobic Unparticle, that would let us evade some of these constraints. So, assuming fi ^ 0, let us first look at 
the typical Primakoff process with a fermiophobic Unparticle that could potentially contribute to supernovae (SN) 
cooling 

7(fci) + 7(fc2) ^ Oip) 

Here O is a fermiophobic scalar Unparticle and one of the photons is sometimes assumed to be off-shell. We are 
only interested in estimating a lower bound for the scale breaking parameter fi that would not violate supernovae 
constraints. Thus for simplicity we assume that the photons are transverse and ignore any plasmon effects (where the 
photon gets a longitudinal polarization in the plasma). Using 

2\ /d2 ,,2\J-2 
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and the Feynman rule from Eq. (j4]) we may compute the invariant amplitude to get 



where s = (^1 + ^2)^ = and O is the Heaviside function. From the above amplitude we may estimate the Unparticle 
Primakoff cross-section to be 



Using the above expression, the emissivity may be calculated from 

£^5^0(7 + 7^ CI) 



{nl^ 0-^7^0(7 + 7 ^U) Ecm) 



PSN 

and gives 

d^ki 2 f <Pk2 



■ , 1 [ d-^ki 2 f 



2 



(27r) 



(42) 



2£;i£;2 8A^^' 

where is the average supernovae core density and Tsn is the mean core temperature in the supernovae. We define 
in the usual way the dimensionless variables xi = Ei/kT<^^ and X2 = ii'2/fcTsN to write the emissivity as 

In the limit of perfect scale invariance, — 0, we recover from the above the expressions in may 
be compared with the fermiophilic case of a vector Unparticle ilAv), with fj, ^ 0, coupling to fermions. In this case a 
process such as w — > Uy in the supernovae core can lead to cooling and the emissivity is given by [l6j 

• f ^ ^jj ^ alA.Ti^' p (4xiX2p(xi+x2) fi'y 



Note that the power law dependence of Tsn on the scaling dimension j is very different, between the fermiophobic and 
fermiophilic cases, even if one assumes that the other parameters /i ~ /iy, ~ A^ and 0(1) coefficients c ~ gv But 
we shall see below that the bound on the scale breaking parameter comes out to be of the same order of magnitude 
in both the cases due to the dominance of the Boltzmann suppression factor. 

To be consistent with supernovae models and constraints from SN1987A we require [i^l 

esN < 10^' J/Kg. s 

Eq. (|43p may be integrated and expressed in a compact form in the limit of 3> with the Boltzmann suppression 
term factored out. An explicit numerical computation with the standard values 

Ts„ ~ 30MeV ; psN ^ IQi*' Kg/m^ 

and the above emissivity criterion gives 

^l > 1.25 GeV (44) 

This lower bound for the effective mass, of a fermiophobic scalar/pseudo-scalar Unparticle coupling only to photons, 
is found to be very close to the result in the fermiophilic case ly + ^ Uv where it was found that [IB] 

/ly > 1 GeV 

Based on the bound in Eq. pij) we may take a minimal value for the effective Unparticle mass to be /i ~ 2 GeV. 
This choice already leads to an emissivity well below the allowed limit. A crude upper bound on the scale breaking 
parameter is given by the requirement that p, <^ An. Since we have adopted the prejudice that the Unparticle 
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scale is probably close to the electroweak-symmetry breaking scale this implies that in our model we are assuming 
/i< 0(1) TeV. 

Let us now turn our attention to some of the constraints from cosmology. For the fermiophobic Unparticle sector 
not to interfere with BBN we must require that the Unparticles stay decoupled during that epoch 45]. In the case 
of a fermiophobic Unparticle with ^ — > the condition Fsm-s-o < H (during the radiation dominated era) , by simple 
dimensional analysis, becomes 

^'^'^ A^J' ^10i«GeV, 

This implies that in the range j € [1, 2) we are considering the rate V-^^q red-shifts faster than the Hubble parameter 
[H). It is required that the fermiophobic Unparticles decouple before BBN and not get reheated during SM phase 
transitions to satisfy pu ^ Psm- This may be achieved by requiring that decoupling happen before the QCD phase 
transition at an energy T > 1 GeV [1^. Note that a fermiophobic Unparticle sector, coupling to photons, does 
not generally re-couple after BBN with SM fields because T^^q red-shifts faster than the Hubble parameter. This 
may again be contrasted with the fermiophilic case of a vector Unparticle coupling to fermions through the effective 
operator 
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where re-coupling is possible when 1 < j < 3/2 |45| . 

When scale invariance is broken (/i 7^ 0), since To ~ "-eq (fbl): the relevant processes are Boltzmann suppressed 
by factors of e~^l'^ when fj, > T. Thus the BBN constraints can be evaded as long as /i is above Aqcd ilSl . We must 
point out that we have tacitly assumed that the fermiophobic Unparticles can decay into SM fields [48| when /i is 
sufficiently non-zero and are therefore not stable. There has been some discussion in the literature on this particular 
issue [8, 46, 48]. All the arguments above suggest that while calculating the Uehling shifts with p ^ we must also 
consider the possibility that fi is higher than the minimal value of 2 GeV. Thus without loss of generality we will 
consider a range 

2 GeV < fi <t: Mz 91 GeV 

for both the scalar and pseudo-scalar Unparticles such that there is still a substantial conformal window (/x, Ajj]- The 
above choice would also ensure that any modification to the gauge kinetic term (Aa~^) near the scale /i is within 
experimental limits [l2| and that the effects due to SM Higgs- Unparticle mixing, if present, are suppressed p^. 
We will demonstrate in a short while that the actual energy shifts are relatively insensitive to small shifts in the fj, 
parameter. 

The pseudo-scalar Unparticle potential energies for two choices of the scale breaking parameter /i are shown in 
Fig. [5] This figure is to be compared with that in the case of perfect scale invariance illustrated in Fig. [3] It is noted 
that the potential in the /i ^ case has been suppressed by 0(1) factors compared to the /i — > case. The general 
features of the Unparticle Uehling potential nevertheless remain unchanged between the /i — > and /i ^ scenarios. 

Let us now calculate the energy shifts in the muonic lead transitions when /i 7^ 0. Consider the pseudo-scalar 
Unparticle with scale invariance broken. We may re-write Eq. (I40[) again by defining the dimensionless variable 
z = |q| r. The expression becomes after simplification 



647r4A^-' sm(j7r) J Jo Jr 







1 + log 1 + } + log 



^2 j,2 y 1^ ^2 ^2 J ^2 ^2 1^ J^J^2 



(1-X)/X2 



where we have made an over-subtraction at q = to get a consistent q — >■ limit. From the above expression we 
expect that for any fixed value of the radial coordinate (r) the dominant contribution to the z-integral should come 
from the region of integration with 

o„2 

< (46) 

or its vicinity. For the muonic lead system we are primarily interested in a suitable value of the nuclear charge density 
parameter Q, in Eq. (|23p . was found to be C ~ 4fm^. We found that a minimal value of /i satisfying astrophysical 
constraints is /i ~ 2 GeV = 10.14 fm^^. From these observations we note therefore that for a fixed value of r 
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FIG. 5. The Unparticle Uehling potential energies for the pseudo-scalar case with the same parameters as Fig.|3]except assuming 
is 2 GeV (left panel) and 10 GeV (right panel). The j — 1.01 (dashed line) and 1.15 (small-dashed line) potential energies 
are again illustrated as functions of r. Note that the magnitude of the potential at any specific value of r has been reduced in 
both cases compared to the case of perfect scale invariance ^ — > 0. 



This means that in the relevant region of the parameter space we may expand the logarithm in Eq. (j45p as 



log 1 



1 



.2 N 2 



^2 ^2 J ^2 J.2 2 y^2^2 

to get an integral that is separable in x and z variables. Note that this approximation becomes more and more 
accurate as we raise the value of /x. The variable separable integral may now be evaluated analytically or numerically 
to calculate the expression in Eq. (|45|) . The calculated energy shifts to the 15 — 2P transition for various values of 
the scaling dimension j and /i are shown in the left panel of Fig. [51 It is noted that compared to the — >■ case 
the energy shifts are smaller, assuming all other parameters remain the same, by 0(1) factors. Thus we find that 
incorporating broken scale invariance with a non-zero value of /i does not seem to alter the energy shifts drastically 
from their /i ^ values and the changes are only by whole number factors. Once again, maintaining continuity, the 
ordinary pseudo-scalar case is recovered as j approaches 1 due to Eq. (|20p . If /i is very large and the corresponding 
energy shifts very small then the finite contributions from the higher order counter terms in the effective Lagrangian 
may become important and the numerical approximation we adopt, of keeping only the lowest order terms from ^ 
and ([5|), may break down. 

In the scalar Unparticle case the expression for the energy shift becomes 

^\^,-^-'^r,,^^,-^'''^''> (47, 

G'ln'^A-/ sm(j7r) J Jq ] r Jq z 

after an over-subtraction at (7 = 0. In the relevant region where x z"^ < 2r^/C ^ /i^ we may again simplify the 
integral and compute it to obtain the energy shifts. The calculated values are shown in Fig.|6l right panel, for various 
values of the /i parameter. The limit of j — >■ 1"*", as before, corresponds to the ordinary scalar case. Again, the 
correction from the scalar Unparticle vacuum polarization is found to be larger in magnitude than the corresponding 
pseudo-scalar Unparticle case. But this may as before be an artifact of our approximation in Eq. (|20p . The corrections 
are both positive as in the /i — case. 

Once again we observe that the variation of the scalar Unparticle Uehling shift with is only by C(l) factors. 
This implies that the energy shift is relatively insensitive to changes of the scale breaking parameter, across a wide 
range, in both the scalar and pseudo-scalar cases. More specifically the Uehling shifts are relatively unchanged even 
for /i ^ m^- ~ 0.1 GeV. In this sense the Unparticle oblique correction is sensitive to a wide range of ^ values within 
our approximations. Let us try to understand this a little better. 

Let us look at the expression 
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FIG. 6. Magnitude estimates of the energy shift in the IS — 2P muonic lead transition due to pseudo-scalar (left panel) and 
scalar (right panel) Unparticle vacuum polarizations. The plots are for ^ ~ 2 GeV, 6 GeV and 10 GeV. With these choices of 
^ there is still a substantial conformal window between /i and Ku- The other parameters are taken to be the same as in Fig. 2] 
Note that the energy shift is lower in both cases compared to the /i ^ case, but not drastically. In fact it is observed that 
the energy shifts are relatively substantial even for /i S> rn^- ~ 0.1 GeV compared to ^ — >■ 0. We will explore the reasons for 
this insensitivity to ^ shortly. In general it is observed that as we increase the effective mass the energy shift decreases. As 
before the energy shifts should be interpreted as accurate only up to undetermined 0(1) factors. 



from Eq. (j45|) for the pseudo-scalar Unparticle that we have re-written to include all the dependences on ^. The first 
term in the above expression tends to zero in the /i — >■ limit and is absent in the case of perfect scale invariance. 
Even in the /i ^ limit we note that since x 2r'^/Q <^ fi^ in the relevant region the first term is relatively 
suppressed. Thus the correction from the first term when /i ^ is small, even though the corresponding term in 
the case of — > is completely absent. Now let us look at the second term in the above expression, specifically the 
/i^/M^ factor inside the logarithm. Again that factor inside the logarithm is obviously absent when — ?> 0. Since 
we have adopted the ansatz that the Unparticle scale is in the vicinity of the electroweak scale we haveup to 0(1) 
factors M ~ Ay ^ v. Thus at the level of our approximation for typical values of allowed we have /x^/M^ ^ 1. 
This again implies that the correction in the /i 7^ case, with respect to the — > case, from the second term is not 
very drastic. 

For the scalar Unparticle case the arguments proceed exactly as above for the relevant expression 
Jo z 

and once again we conclude that the corrections to the /i — ^ case from the additional factors are not very large. Thus 
the relative stability of the Uehling energy shifts to variations in the scale breaking parameter /i may be traced to the 
momentum cut-off imposed by e~^^ leading to the observation in Eq. and the fact that /i < M ~ Ku ^ v 

for typical values. 

The typical values for the Uehling energy shift in Fig.|6]are in the range 0(0.1) cV - 0(0.01) cV. The O(0.1)eV 
as we commented previously is comparable to the contribution from light-by-light scattering in QED. One may get a 
feel for the lower value 0(0.01) eV in the range by noting that it is of the same order of magnitude as the corrections 
from the QED fourth-order Lamb shift, to higher angular momentum states in /i^ — Pbg2^ [s^. For higher angular 
momentum states the muon's anomalous magnetic moment induces an additional spin-orbit interaction in the muonic 
atom and this contributes to an energy shift. For instance in the /i~ — Vh^^ bg — Af transition the leading contribution 
to the fourth-order Lamb shift at order (Za) was estimated to be |34| 

AE^-Jf^j ~ 0.025 eV 

The values of the Uehling shifts in Fig. [5] may again be compared to the uncertainties in the precision measurements 
of 2pi/2 — lsi/2 and 2^3/2 — lsi/2 transitions in Table IV and the discrepancy between theory and measurement for 
these transitions in (|35p . It is noted that for typical values of the model parameters {j, A^, c and ^) the Uehling shift 
is again about a factor of lO'^ — 10^ below the values in ([55)1 . 

We may also calculate the Uehling shifts for the low-lying / = 1 and I — 2 states with respect to IS". It is observed 
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that for /i 7^ the Unparticle scalar and pseudo-scalar corrections again follow a hierarchy 



consistent with our expectations in Eq. (|33|) . for the case of perfect scale invariance, and also with our computations 
for the /X — case before. We conclude that the choice of /x 7^ does not change the level structure of the Uehling 
shifts. Moreover the variations are generally of the same order of magnitude as in the /i — > case. For example, 
assuming /i ~ 2 GeV and Ajj ~ v, the variations between the I ~ and 1 states for n = 2 are again of the order of a 
few meV and for n = 3 of the order of a few 0.01 meV. 

A speculation is that one may do precision spectroscopy of the low-lyin g at omic states in parallel with proposed 
efforts to observe coherent muon- electron conversion in muonic atoms (see |49| and references therein). The relevant 
process in coherent muon-electron conversion is 

+ N ^ + N 

where iV is a nucleon. It is believed that probes of resonant muon-electron conversion near a nucleus may be able to 
achieve a higher sensitivity to lepton-flavor-violation (LFV) compared to direct conversions jsoj ] 

^ + 7 

This opens the possibility that one may also perform precision spectroscopy on low-lying muonic atom states in 
these forthcoming experiments and reduce some of the discrepancies in Eq. (|35|) . But reducing the discrepancy to a 
level of 0(0.01) eV looks very improbable to us. 

Since muon-conversion is a coherent process one might expect that the probability of muon-conversion in an atom 
would go like ~ (or ^ A^). So when normalized to the muon-capture cross section we have heuristically 

^ T[fi- + {A,Z)^e- + {A,Z)] Z^^ ^ 

r [fi- + {A, Z)^v^ + (A, Z-l)]^ Z ^ 

while for the Uehling shift we are interested in, as we noted previously, the dependence on the atomic number Z goes 
hke 



It is not clear that the muonic lead system, with Z = 82, that we are considering is suitable for the LFV measure- 
ments since it has been shown through detailed calculations [5l| that for coherent muon-electron conversion the most 
ideal range is Z € [30, 60]. Thus it would be interesting, as we have previously mentioned, to explore the possibility 
of measuring the scalar/pseudo-scalar Uehling shifts in other atomic systems, with an intermediate Z value, where 
the energy shifts may still be substantial while the system is also of interest to coherent muon-electron conversion 
experiments. It would also be interesting to estimate the Unparticle vacuum polarization effects in muonic atoms 
where high-precision LASER spectroscopy might be possible [52|. We hope that this work will be a modest pointer 
in this direction. 

It is of course an open possibility that if for some reason the interaction energy scales (A-y, A-y) are too large (^ 
0(1) TeV) or the coefficients (&, c) in Eqs. (|4]) and (O are very small the energy shifts will be even more suppressed and 
the Unparticle Uehling shift, even if it exists, will be nearly impossible to detect. Even under optimistic assumptions 
about the model parameters it is possible that the theoretical difficulties in calculating the required higher-order 
QED/nuclear effects in muonic atoms may be insurmountable or that obtaining precision spectra of the low- lying 
states is very difficult. In such a scenario the only hope would be to look for fermiophobic Unparticles in very high 
energy colliders or other systems that are more sensitive to fermiophobic Unparticles. 

IV. SUMMARY 

In this work we tried to study some of the probable effects of a fermiophobic scalar/pseudo-scalar sector on bound 
state energy levels, specifically low-lying muonic atom levels, as a consequence of oblique corrections to the photon 
propagator. 

Considering the scalar and pseudo-scalar fields to be Unparticle operators, without loss of generality, we examined 
the functional forms of the vacuum polarization functions and the induced Uehling potentials. Some interesting 
theoretical observations were made on the singular nature of the Unparticle induced Uehling potential and the behavior 
of the energy shifts in the limit of the scaling dimension approaching unity. 
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It was estimated that for an Unparticle scale near the scale of electroweak symmetry breaking, in the low TeV range, 
the energy shifts in the low-lying muonic lead transitions could typically be of the order of a few 0.1 eV to a few 0.01 
eV for some natural values of the model parameters. It was also pointed out that these magnitudes are comparable 
to bound state QED corrections, to the higher orbital angular momentum transitions in muonic lead, from the virtual 
Delbriick effect (light-by-light scattering) and the fourth order Lamb-shift (at order a^(Za)) respectively. These 
conclusions are relatively unchanged even when one incorporates a breaking of the scale invariance by introducing an 
effective Unparticle mass fi. 

But the current discrepancy between muonic-lead spectroscopy and theory, especially nuclear theory, makes an 
interpretation of the Unparticle Uehling shift, if it really exists, extremely challenging. A conservative estimate is 
that such an interpretation would require an improvement in the discrepancy between theory and experiment, from 
about 20 years back, by a factor of 1000 — 10000. The recent, partial resolution of the long standing discrepancy in 
the A2p and A3p NP calculations with results from muonic lead spectroscopy [i^ |4l| by A. Haga and co-workers [s^ 
is a promising step in this direction. 

We also mentioned that in cases where the UV sector has a very large fermion multiplicity the above contribution 
may be greatly enhanced, but appealing to arguments of naturalness we do not think this to be very plausible. But 
many interesting models being considered today (for example string- inspired QCD-like models ^21i] ) allow for the 
possibility of a large fermionic sector and this perhaps beseeches us not to discard the possibility of fermiophobic 
Unparticle oblique corrections prematurely. 

The other interesting direction is to consider intermediate- Z muonic atoms where the nuclear uncertainties may be 
much better controlled while at the same time have sufficient fermiophobic Unparticle contributions to muonic-atom 
transitions, by virtue of the SE^^"'- / E^^^^^ ^ Z"^ enhancement. This is left for future work. 

In the context of the present study, of a possible fermiophobic Unparticle scalar/pseudo-scalar sector, we also briefly 
considered constraints from astrophysics and cosmology and put bounds on the fermiophobic Unparticle effective 
masses. Finally we speculated on improving muonic lead spectroscopy and theory in the context of forthcoming 
experiments that will study coherent muon-electron conversion. 
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